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Abstract 

We examine fully coherent two-pulse propagation in a lambda-type medium, under two-photon resonance 
conditions and including inhomogeneous broadening. We examine both the effects of short pulse preparation 
and the effects of medium preparation. We contrast cases in which the two pulses have matched envelopes 
or not, and contrast cases in which ground state coherence is present or not. We find that an extended 
interpretation of the Area Theorem for single-pulse self-induced transparency (SIT) is able to unify two- 
pulse propagation scenarios, including some aspects of electromagnetically-induced transparency (EIT) and 
stimulated Raman scattering (SRS). We present numerical solutions of both three-level and adiabatically 
reduced two-level density matrix equations and Maxwell's equations, and show that many features of the 
solutions are quickly interpreted with the aid of analytic solutions that we also provide for restricted cases 
of pulse shapes and preparation of the medium. In the limit of large one-photon detuning, we show that the 
two-level equations commonly used are not reliable for pulse Areas in the 2n range, which allows puzzling 
features of previous numerical work to be understood. 







I. INTRODUCTION 



The study of the fully coherent propagation of light in the quantum domain, where the non- 
classical response of a resonant medium is taken into account, was opened by the famous work of 
McCall and Hahn [1] . They used low-temperature ruby and nanosecond-scale pulses, and reported 
experimental realization of the first optical solitons, in agreement with their novel theory of self- 
induced transparency (SIT). McCall and Hahn supplied the theoretical foundation for the SIT 
effect by self-consistently combining a two-level model for the interacting Cr 3+ impurity ions with 
a semiclassical treatment of the field [2] . The understanding of coherent optical pulse propagation 
in two-level resonant media was rapidly extended Q, U, [f3, |g, Q|. 

Our interest here is in coherent two-pulse propagation (for an early example, see Q]). Many 
different parameters of both the medium and the pulses themselves can play physically distinct 
roles during propagation. We adopt common restrictions made in treating lambda-type media in 
order to provide a close focus on a few interesting questions. For example, we assume that the 
Rotating Wave Approximation is reliable and that the detunings of the two pulses are large and 
are also equal, so that two-photon resonance applies. 

To focus attention on coherent evolution we also assume that the pulse durations are short 
enough to ignore homogeneous relaxation of the atoms, but we allow inhomogeneous broadening. 
To restrict the range of atomic effects we assume that the pulse is too weak to cause excitation 
beyond the three lambda levels, that the "total" pulse Area is of order 2tt, and that the propagation 
path is long enough to allow asymptotic conditions of propagation to be reached. A standard sketch 
is shown in Fig. [TJ Because of our focus on fully coherent propagation there are natural connections 
to SIT, and because our medium is of lambda type, some connections are expected with Harris's 
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y-induced transparency (EIT) [9( and traditional stimulated Raman scattering 



13l | and fully coherent "transient" SRS 14], all of which usually employ a quasi-cu; 



or undepleted pump field instead of a short pump pulse. 

We first present a class of analytic solutions to our two-pulse nonlinear evolution equations 
(three- level Maxwell-Bloch equations) . The solutions support our earlier observation [lj]] that in the 
situations of interest one can divide the evolution process into three widely separated propagation 
zones: a strong interaction or "pulse transfer" zone, and the zones in the distant past and distant 
future of the transfer zone. The analytic solutions are analogs of previously known resonant SRS 
solitons presented by Bol'shov, et al. [16] and our use of the Park-Shin Backlund solution approach 



171 ] allows incorporation of far off-resonant pulses and a medium with either pure or quantum- 
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mixed ground states. We find that our solutions are constrained by a new two-pulse Area condition 
similar to the 2tt Area condition in SIT. 

The analytic solutions are then used in several ways. We have found in numerical modeling 
that the solution formula and the new two-pulse Area condition serve as accurate predictors of 
output pulse shape and total Area, even for a variety of input conditions not conforming to the 
analytic solutions. This is closely reminiscent of the behavior of coherent short pulses in one-photon 
resonant two-level media, which evolve coherently toward the 2tt sech pulse of SIT, even if injected 
into a medium with shape and Area not in that form. 

We also make use of the analytic solution formulas to check departures from adiabatic predic- 
tions. An adiabatic restriction is normal in time-dependent treatments of EIT effects, and we are 
interested in its limits of applicability in the context of propagation. Similar questions can be asked 
about SRS, commonly treated as an adiabatic propagation process governed by approximate far 
off-resonant conditions, where a large product of detuning and pulse duration, rA 3> 1, is used to 
justify an approximate two- level set of evolution equations. We solve numerically for the evolu- 
tion behavior in the large-detuning case, and demonstrate the inadequacy of the reduced-equation 
approach if either pulse has near-unit Area. 

In some early numerical work with two-pulse evolution in three-level media [la ], a striking 
and unexplained pulse breakup was reported for solutions of the reduced two-level equations. 
Our results, obtained in a similar physical situation but with the more fundamental three-level 
equations, suggest a resolution of this effect, as pulse breakup related to higher order SIT- type 
multi-pulse formation. 

Finally, the analytic solutions permit a new kind of comment on two-pulse propagation in 



lambda media. We use them, along wit 



i numerical solutions, to show a strong contrast with 



the earlier reports by Kozlov and Eberly [19] of relative pulse stability of EIT-type and SIT-type 
propagation behavior. We show here that the presence or absence of quantum coherence of the 
ground-level preparation of the medium is the crucial property in determining asymptotic stability, 
and we conclude that mixed-state preparation leads to SIT-type stability, opposite to the conclusion 
reached for pure-state preparation. 

II. THEORETICAL MODEL 

We consider dual-pulse propagation in a medium of three-level atoms in the lambda config- 
uration as shown in Fig. [TJ Given linear polarization, the electric fields of the individual laser 
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pulses can be written as scalars: E a (x,t) = £ a (x,t)e~ l ^ kaX ~ UJat ^ + c.c., where £ a is the slowly 
varying envelope of the electric field, k a is the wavenumber, and u) a is the frequency, with a sim- 
ilar formulation for the second pulse Eb(x,t). For the 1-3 transition we have the Rabi frequency 
£l a = 2d a £ a /h, where d a is the interacting component of the transition dipole moment, and the 
detuning is A = U3 — (u>i + uj a ). We use corresponding notation for the 2-3 transition. Note that 
the same A serves as the one-photon detuning of each transition below resonance, which implies 
exact two-photon resonance. 




FIG. 1: Three-level atom in two-photon resonance with level 1 connected to level 3 via the Rabi frequency fl a , 
referred to as the pump field, and level 2 connected to level 3 via the Rabi frequency Qt, referred to as the Stokes 
field. Each laser field is detuned from the excited state resonance by an equal amount A, and the excited state decay 
rate is given by 73. 

The Hamiltonian of the system in the rotating wave picture is given by 

H = ^A|3)(3|-^|1}(3|-^|2)(3|-^|3)(1|-^|3)(2| 

= RA|3>(3| - (^|1) + ^|2))(3| - |3)(^(1| +^(2|), (1) 

where the factorization in the last line suggests introducing the "bright" and "dark" state repre- 
sentations of the two lower states |1) and |2): 

\ B ) = ?&\l)+}&\2) and \D) = h—\l) - h—\2), (2) 
2 2 17 2 1 1 2 1 " K ' 

where the interaction terms in H clearly depend only on \B), while the orthognal "dark" state \D) 
[20I does not participate in the temporal dynamics. 

Dynamical evolution is confined to \B) and |3) only so long as spatial evolution is ignored. When 
spatial evolution is permitted, and the dark state is a participant, it can be difficult to maintain 
some of the restrictions that are commonly imposed in treating two-field physics (fixed control 
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field intensity, a permanently weak probe field, etc.). Some effects of dark-state participation 
in propagation have been pointed out previously (e.g., see [2l|). and Kozlov and Eberly more 
recently showed that the dark state even completely controls evolution asymptotically 
when propagation occurs in a lambda medium of two coherently prepared ground states (sometimes 
called "phaseonium" [23J]). The contrast extends to stimulated Raman scattering, historically the 
first lambda-medium propagation effect, and not associated with an extended prepared coherence 
between the two ground levels. 

We now derive analytically exact two-pulse solutions for lambda media. This requires an ex- 
tension of previous soliton theories for resonant SRS propagation (see [if], ) • Here we depart 
substantially from resonance, however still near enough that we can neglect anti-Stokes effects or 
interaction with any other excited states. We adopt SRS nomenclature and will refer to £l a as 
the pump field and f2& as the Stokes field. We will focus on non-phaseonium medium preparation, 
but include a variety of mixed-state descriptions of levels 1 and 2. This permits us to account for 
incoherent initial population distributions that may arise from thermal effects, incoherent optical 
pumping of the ground states, etc. The solution method we use allows for significant changes in 
both the pump and Stokes pulses, and gives analytical insight into previously pursued numerical 
studies of short-pulse Raman propagation . 

The von Neumann equation for the atomic density matrix is 

ihp=[H,p], (3) 

where the over-dot signifies d/dt and where H is given in Eq. (pQ). For later reference, the equations 
for the density matrix elements are: 

Pu =i y P3i ~ i-^Pw ( 4a ) 

P22 =*y/032 - i^P23 (4b) 
P33 = ~ l ^P31 + ^P13 - 1 y P32 + i — P23 (4c) 

o o* 

P\2 =iy P32 ~ «y P13 (4d) 

pi 3 =iA/)i 3 - Zy/)i2 +iy(/) 3 3 - Pn) (4e) 

p23 =£Ap 2 3 - iy P21 + (P33 ~ P22)- (4f) 

Without damping terms, Eqs. ([3]) and ^ imply that the duration of the laser pulses under 
consideration is much shorter than all homogeneous relaxation processes (such as the excited state 
decay, with rate 73, shown in Fig. [[]), thus justifying their neglect. 



For simplicity we will change to a travelling reference frame where the time derivatives are with 
respect to the retarded time T = t — x/c and the spatial derivatives are with respect to Z = x/c. 
This specification gives 

d d d d d 

dt = df ^ C d-x + di=dZ- (5) 
The field evolves according to Maxwell's equations. In effect, the interaction Hamiltonian 
itself, which is just the dipole interaction d ■ E(x,t), evolves according to Maxwell's equations. 
We assume that the pump field only interacts with the 1 — > 3 transition, the Stokes field only 
interacts with the 2 — ► 3 transition, and the fields move only in the positive x direction. These 
assumptions are justified by the inequality \{u2 — u>i)/A| 3> 1 together with the slowly varying 
envelope approximation, and they give two distinct field equations from Maxwell's equations: 

= ~^aPl3 (6a) 

— = ~lfi b P23- (6b) 

Here N is the density of atoms, fi a = Nd a u>^/heo is the atom field coupling parameter for the 1 — > 3 
transition, and similarly for transition 2 — ► 3. For simplicity in the derivations to be presented, 
we have not included inhomogeneous broadening in these equations, but it will be incorporated 
in the end. We note that when fj, a = ^ = fx these results can be combined into a single matrix 
equation for the spatial evolution of the Hamiltonian: 

9Z=-iiW,p], (7) 

where 

/l \ 



w = - 1 - 

2 



I) 1 , (8) 

v° ° -v 

and this matrix format will be useful shortly. For notational simplicity we have suppressed the 
functional dependences of the density matrix elements in Eq. ([3|) and the Rabi frequencies in Eqns. 
J6} on T and Z (or equivalently x and t), and we will do so for the remainder of the paper. 

III. PARK-SHIN BACKLUND SOLUTION METHOD 

Eqns. Q and Q are known as the Maxwell-Bloch equations for the three-level lambda system. 
For comparison with adiabatically approximated solution methods as well as with numerical meth- 
ods it is highly useful to have solutions both in analytic form and without adiabatic assumptions. 
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The inverse scattering technique has previously been used to construct resonant solutions 
These equations are integrable when fJ> a = l^b = Hi which we will assume from now on. We will 
solve these equations using the Park-Shin (PS) formulation [l^] of the Backhand solution method. 
We will briefly sketch the PS method, and slightly generalize their approach to include arbitrary 
detuning of fields. 

One begins by finding a pair of linear operators L\ and L 2 , which will allow for a linear represen- 
tation of the nonlinear Maxwell-Bloch equations when subject to the constraint that [Xi,!^] = 0. 
We define these operators as: 

Li = A + U + XW (9a) 
L 2 = (X-A)-^ + V, (9b) 

where A is the spectral parameter, and U and V are given by 

U = ±H° (10) 
n 

v = (11) 

where H° is the Hamiltonian in Eq. ([I]) with A = 0. To recover the Maxwell-Bloch equations one 
imposes the commutator relation [Li,L 2 ] = and then compares terms with like orders of A. In 
doing so one finds the two independent 3x3 matrix equations 

ih^ = [H,p] (12a) 

% = \ W ^ (12b) 

Eq. (|12ap is simply the von Neumann equation, and it is straightforward to show that Eq. (|12b[) 
is equivalent to the matrix form of Maxwell's equation presented in Eq. ([7]). Thus we see that 
the commutator [Li,L2] = correctly reproduces the Maxwell-Bloch equations, and is equivalent 
to the integrability condition, % - % + [V, U'] = 0, of the "Lax" Pair U' = -(U + XW) and 

The PS method now proceeds as follows. Because L\ and L2 commute, they have common 
eigenvectors, which we represent here in the form of matrix M: 

L\M = (13a) 
L 2 M = 0. (13b) 
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We now apply the Backlund transformation to Eqns. (|13p by writing the solution matrix M in 
terms of a known "proto-solution" matrix, which we denote as M^°\ and a dressing operator 
<3? : M = $M(°). The proto-solution can even be a trivial solution. We will write $ as 

$ = 1 + 4(2P-1), (14) 
A 

where P is a projection operator in the 3x3 atomic space, implying the existence of a vector \s) 
such that 

„_I«X*I (15) 



(s\s) ' 

and T] is a parameter of the Backlund transformation, which we take to be real. Eqns. (113p and 
(|14p combine to give: 

d 

— $ - <1>(Lt( ) + \W) + (U + \W)<S> = (16a) 
(A- A)— <5> + V<S>-$V i0) =0, (16b) 

where l/( ) and F< ' are the known solutions of the fields and density matrix making up the proto- 
solution matrix M^°\ We now insert the definition of into Eqns. (|16p . and compare terms with 
like orders of A, which results in four independent equations. This allows one to eliminate equations 
involving the space and time derivatives to arrive at the following two equations: 

U = If® -2irj[W,P] (17a) 

v -^{ 2p - i -'^ v '°'( 2p - 1 + ^)- (17b) 

Eqns. (I17p relate the unknown solutions given by U and V to the known solutions given by E/(°) 
and and the matrix P. 

The only unknown quantity in Eqns. ()17l) is the projection operator P (or the vector \s)). 
We can find an equation for P in terms of known quantities by again using Eqns. (I16p . but this 
time solving for terms containing the derivatives. The result, after simplifications made possible 
by utilizing projection operator properties, is two independent linear differential equations for the 
vector \s), given by: 

(Jf + C/(0) ~^ W ) I s ) =0 ( 18a ) 
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Since the commutator [U^ — irjW, V^] vanishes, the solutions to Eqns. (|18ap and (|18bj) can be 
exponentiated together to give 



exp (-U^T + ir/WT + ^—^V^Zj \u), (19) 



where \u) is an arbitrary constant vector, with coefficients (ui,U2,U3), which is specified by the 
initial condition of the pulse. At this point, one can introduce inhomogeneous broadening by simply 
Doppler averaging the third term in Eq. (|19|) with the Doppler distribution function 

F(A) = JpLe-CA-Am*^ (20) 



where is the inhomogeneous lifetime and A is the line-center single-photon detuning of both 
laser fields. To be consistent one must also average the right hand side of Maxwell's equation (|6|) 
with the same function. This could have been incorporated from the start. 

Finally, since we are interested in the actual density matrix p and the fields £l a and (instead 
of U and V), we write them in terms of the vector \s) using Eqns. (|17b|) and (|17ap . The density 
matrix solution is: 

rj 2 f n \s)(s\ V + iA\^ 0) A >)(s| rj-iA \ ^ 



rf + A 2 V, (s\s) r\ ) r \ (s\s) 7] 

and the individual field solutions are: 

= 2i(l|l/|3) = Ol°) - 4 ? r,^M (22a) 

(s|s) 

Q b = 2i(2\U\3) = UP - 4n / 2| ^,y (22b) 

(s|s) 

To summarize, the coupled non-linear Maxwell-Bloch equations (J3]) and © have been re-written as 
two coupled linear differential equations (|18ap and (j!8bj) and their solution is given in equation (|19p . 
The solution depends on a known proto-solution to the original coupled Maxwell-Bloch equations 
given by and V^°\ or equivalently Q^, , and p(°\ Any solution of the Maxwell-Bloch 
equations can be used as this proto-solution. In the next section we will show how even trivial 
proto-solutions will yield interesting results. 



IV. TWO-PULSE ANALYTIC SOLUTIONS 

We now use the PS approach to solve the MB equations explicitly for the case off-resonant SRS 
pulse transfer. The solution method allows one to treat mixed as well as pure-state media and 
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also allows for treatment of spatially-dependent initial level populations 24j, |25j. The choice of 
proto-solution gives the initial condition of the medium and pulses at T = —oo. 

The state of a lambda atom with no population in the upper level has the form |^ a tom(0)) = 
ci 1 1) +C2 12) , so if we write c\ = ae 1 ^ 1 , and C2 = /3e"^ 2 , with real a and j3 constrained by a 2 +f3 2 = 1, 
we can identify phase effects easily in the corresponding density matrix 





p(0) = |V>atom(0))(Vatom(0)| 



a 

af3e~ 




(23) 



af3e ie 

o oy 

where 5(f) = c/>i — <j>2- Complete initial dephasing is a natural situation consistent, for example, with 
thermal mixing between the two ground state populations or with incoherent optical pumping 
between them. In this case the off-diagonal terms of the density matrix vanish and we have a 
diagonal initial density matrix: 

(a 2 (A 



P(0) 



V 



(24) 



p 2 
0/ 

We easily see that the state specification given in Eq. (|24j) . along with the zero values = 



and 



(0) 



i. — 0, constitutes a trivial but exact solution of the MB equations. It is also useful because 
the PS Backlund approach constructs a non-trivial solution from it. The trivial solutions give 
CA°) = and = —igp(0). Since these solutions give a diagonal matrix in the exponent of Eq. 
(|19p . the components of vector |s) are easily calculated. They are: 

n/„ o ,/ i 

2r . 
1 

2t. 

T 
2t 



(l|s) = exp 
(2|s) = exp 

(3| S ) : 



2 2 

a [it 



a 2 ^ 



1 + iAr 
1 

1 + iAr 



zexp 



(25a) 
(25b) 
(25c) 



where we have identified 1/rj = t as the nominal pulse width and take u\ = 1, U2 = 1, and U3 = — i 
for greatest simplicity. Eqns. ([22]) and ([25]) combine to give the two pulse solutions: 

ia 2 5Z 



a, 



Ae 



4e" 



2cosh(T/r - o> 2 kZ) + exp (T/r + (a 2 - 2/3 2 )kZ) 
2cosh(T/r - (3 2 kZ) + exp (T/t + (/3 2 - 2« 2 )kZ) 



where the functions k and 5 are given by the inhomogeneous average: 

2 J A-i/r 



(26a) 
(26b) 

(27) 



or 



F(A)dA 



and 5 = — 



AF(A)dA 



(28) 



2r7_ ooA 2+(I)^ 2J_ x A2 + (I) 1 

Here 5/c adds to the index of refraction and k/c is the absorption depth for weak-field excitation 
and so it sets the distance scale for propagation effects. 

Substantial transient atomic coherence develops during the course of propagation and off- 
diagonal density matrix elements become non-zero. These solutions are calculated from Eqns. 
(J2H) and ([25D, giving: 



1 



Pll 

P22 
£»33 
P\2 
Pl3 
P23 



1 + (At) 2 
1 

1 + (At) 2 
1 

1 + (At) 2 
1 

1 + (At) 2 
1 

1 + (At) 2 
1 



a 2 (|/n| 2 + (AT) 2 )+/3 2 |/i 2 | 2 
(a 2 |/i 2 | 2 + /3 2 (|/ 22 | 2 + (AT) 2 ) 

(a 2 |/i 3 | 2 + /3 2 |/2 3 | : 
(a 2 (f u -iAr)f 12 +p 2 (f 22 + iAr)fr2 
(a 2 (fn-iAT)f 13 + p 2 f 12 f 23 : 

(« 2 A* 2 /i3 + /3 2 (/22-iAT)/ 23 ) 



1 + (At) 2 

where the space-time dependences are contained in the functions fif. 



hi 

/22 

/12 
/l3 
/23 



2 sinh(T/T - a 2 nZ) - exp (T/t + (a 2 - 2/3 2 )kZ) 
— 2 cosh (T/t — a kZ) + exp (T/t + (a 2 - 2/3 2 )acZ) 

2e -i(^-^)5X e T/r-^^ /jD(Z)T) 



T>(Z, T) 
D{Z,T) 



lie 



2 8Z 



/D(Z,T) 



2ie ~^ 2sz e^ 2 ~^ z /D(Z,T), 



(29a) 
(29b) 
(29c) 
(29d) 
(29e) 
(29f) 



(30a) 
(30b) 

(30c) 
(30d) 
(30e) 



and the denominator function D(Z, T) is given by 



D(Z,T) = 2 cosh (T/t - a 2 nZ) + exp (T/r + (a 2 - 2/3 2 )kZ). 



(31a) 



A. Input Regime: -kZ ^> 1 

The general pulse and density matrix solutions given in Eqns. (I26p and f)29[) are cumbersome 
and resist rapid interpretation. However their key asymptotic features are clear. First, in the far 
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negative spatial region, i.e. — kZ 3> 1. The asymptotic forms of Eqns. (|26[) simplify to 

2g— ia 2 8Z 

n a = sech(r/r - a 2 nZ) (32a) 

r 

n b = 0, (32b) 

which one can easily recognize as the McCall-Hahn 2ir hyperbolic secant pulse solutions [l| for SIT 
on the 1 — > 3 transition, and no pulse on the 2 — ► 3 transition. We interpret this as an "input" 
pulse advancing with a group velocity given by Vg Jc = (l+a 2 Acr) — 1 . The density matrix solutions 
also simplify substantially in the same limit. They are 



p n = ^_^ a 2 (tanh 2 (T/r - a 2 nZ) + (At) 2 ) (33a) 

1 + (At)^ ' 

P22 = P 2 (33b) 

P33 = i , 77 n 9 » 2 sech 2 (T/r - c?kZ) (33c) 
1 + (At)' 2 

Pl2 = (33d) 

Pl3 = i T+ (Ar) 2 " 26 "" 2 ^ sech(T/T - a 2 «;Z) ( tanh(T/V - a 2 KZ) - i At) (33e) 

P23 = 0. (33f) 



One can see that all of the interaction in this asymptotic limit is localized to the 1 — * 3 transition, 
and the three-level solutions simplify to those for a two-level atom. 



B. Output Regime: nZ ^> 1 

Similarly we also look at the solutions for large positive propagation distance where kZ 3> 1, 
the "output" regime. These solutions are given by 

n a = (34a) 

2e -ii3 2 sz 

n b = sech(T/T - (3 2 kZ), (34b) 

T 

which is another SIT solution except that the excitation has been transferred to the Stokes pulse, 
and here fl b moves with group velocity Vg jc = (1+/3 2 «;t) _1 . The density matrix solutions similarly 
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change. They are now given by 

pn = a 2 (35a) 
P22 = 1 2 /? 2 (tanh 2 (r/r - (3 2 kZ) + (At) 2 ) (35b) 

P33 = YTJKtY ^ S6Ch2 (T/T " ^ KZ) (35C) 

pi2 = (35d) 

Pis = (35e) 

p 23 = i 1 + (A T )2 ^ 2e ~ i/32gZ sech ( T / T - /5 2kZ ) ( tanh(T/r - /3 2 kZ) - zAr) . (35f) 

C. Pure State Input/Output Regimes 

As a limiting case we take the initial preparation of the atoms to be in a pure-state by setting 
a = 1 and (3 = 0. In this case the input pulses for — kZ » 1 are given 

n a = sech(T/r - kZ) (36a) 

r 

n b = 0, (36b) 
while the corresponding input regime excited state probability is 



P33= l + (Ar)2 Sech \-- kZ I- (:57) 



The output pulses for kZ 3> 1 are similarly found and given by 



n a = (38a) 

n h = - sech(T/r), (38b) 
r 

and the output regime the excited state probability is simply 

p 33 = 0. (39) 

One can verify these forms by comparing them to the mixed-state case, taking a = 1 and (3 = 0, 
and applying the appropriate limits for the input and output regimes. 

In this limiting case the input pump pulse is identical in form to an SIT pulse and is coupled 
to the 1 — ► 3 absorbing transition. It causes coherent atom excitation, and travels with reduced 
group velocity Vg^ /c = (1 + kt) _1 . Meanwhile the output Stokes pulse, while similar in temporal 
shape to an SIT pulse is completely decoupled from the medium, and thus moves with group 
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velocity v g /c = 1. This is due to the fact that we have neglected any anti-Stokes interaction, and 
have limited our atomic system to three levels. This is a trivial example of a Dark State where 
lack of population in state 2, rather than quantum interference effects, causes the medium to be 
transparent to the Stokes pulse. 



V. PULSE AREA AND TRANSFER ZONE BEHAVIOR 



Since the solutions in both forward and backward asymptotic regimes reduce to two-level soli- 
tons, in these regimes the pulses are expected to obey the two-level Area Theorem Qj. The Area 
of a pulse is defined to be 

/oo 
n(Z,T)dT, (40) 
-oo 

from which the pulse Areas of solutions (|26p can be shown to be 

e " iz) = W) md 6 » [z) = Wv (41) 

where 



h(Z) = Vl + e 2 ^ 2 -/ 32 )^. (42) 

The depth of medium required for significant change in pulse Area obviously depends on the so- 
called "Raman inversion" a 2 — (3 2 , and this is indcated in the plots of 9b(Z) in Fig. [2j 

From the analytic formulas we derived one easily finds the =fZ asymptotic limits: 9 a (Z) — > 2ir 
and 9b(Z) — ► 2tt, as noted in the previous Section. However, in addition, one can extract a 
surprising new result if one defines a "total" Area 0(Z) for any position Z: 



9(Z) = <j9 2 (Z) + 9 2 (Z) = 2ir. (43) 

This makes the connection to SIT very close indeed, and shows unexpected continuity of total Area 
through the non-SIT transfer zone. The significance of this exact result for us is that the pulse 
strengths and pulse durations of the exact solutions are fully linked through the two-pulse Eqn. 
(|43|) , This gives the impression that the analytic solutions apply only to situations in which at 
least one pulse is "strong", with an Area greater than ir. We will comment further on this point 
in the next Section. 

Another point is that our solutions describe a single-peak pump pulse and a single-peak Stokes 
pulse. For larger total pulse Areas (e.g., 3ir < 9 < 5ir) one would expect to obtain double- 
pulse soliton solutions. We do not consider higher than single soliton solutions here, but they 
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FIG. 2: Curves showing the way the analytic formula (j41j) predicts 0f, to grow toward Z = where its Area 
matches the pump pulse Area 9 a - The distance required for this growth from the initial value 0b — 0.0057T 
depends on the so-called Raman inversion a 2 — 1 , and three examples are given as indicated to the right 
of each graph. These distances are designated as Zt in the text, to indicate the distance of propagation in 
the medium required to reach the "transfer" point where 6b — 6 a . The large plot shows the early growth 
stages for very small values of 9b while the small plot shows the entire growth of the Stokes pulse. 



can be obtained by using a nonlinear superposition rule that can be derived from the Backlund 
transformation. For more we refer the reader to 1171 1. 

n 

We have already reported [155] that the analytic solutions describe the transfer of a 2ir Area 
sech-shaped pulse on the 1 — ► 3 transition to a 2tt Area sech-shaped pulse on the 2 — > 3 transition, 
through an intermediate interaction regime or "transfer zone". We can obviously use this three- 
regime language here, where regimes I and III contain the asymptotic input and output pulses 
described by Eqns. (|32ft - (|35p . and regime II contains the transfer between them, where the 
full three level equations (|26p and (|29p are needed. Regime II can be thought of as exhibiting a 
strong-field SRS process. 

The analytic solutions for an off-resonant interaction (i.e. At = 10) are shown in Fig. [3l and 
the three-regime behavior just discussed is clearly evident. The pulses are plotted as a function of 
the laboratory frame position variable x in units of k/c, and each frame corresponds to a time t 
in units of r. As stated previously these are related to the coordinates presented in the solution 
by T = t — x/c and Z = x/c. In these examples, during the input stage of propagation, only the 
pump pulse is significantly present, corresponding to regime I (frames 1 and 2 of Fig. [3]). During 
"transfer" from pump to Stokes both pulses are intense, indicating that we are in regime II (frames 
3-5). Finally the asymptotic output emerges and only the Stokes pulse is significantly present 
(frame 6) indicating we have reached regime III. 
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amplification. The input pulse on the 1-3 transition amplifies a weak probe pulse on the 2-3 transition. 
Parameters: a 2 = 1.0, (3 2 = 0.0, r w 3T 2 *, and At = 10. 



VI. FLUX CONSERVATION 



It may appear "obvious" that energy conservation cannot be satisfied through the Raman 
exchange regime II. If the original input pulse is composed entirely of type "a" photons as seen in 
Eqns. (|36|) . while the output pulse is composed of different frequency type "b" photons as shown 
in Eqns. (138 j) . the energy difference associated with the different pulse frequencies appears to be 
lost, since the medium stores no energy at the end. 

In the situation at hand the interacting-system conservation law is one of flux, not of energy. 
An expression for Poynting's Theorem in this one-dimensional example can be derived from the 
coupled Maxwell-Bloch equations and © and is given by 



A ( | fia | 2 + |^| 2) + 2 ^^3) =0 . (44) 
In the input regime I where only a pump pulse is present, Poynting's Theorem becomes 
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and for the output regime III where only the Stokes pulse is present we get 

^l^| 2 = 0, (46) 

where the ^33 term also vanishes as shown in Eq. (j39[) since we are considering the pure-state for 
this example. 

The regime III Stokes pulse solution, given in Eq. (|38p . trivially satisfies Eq. (|46p . since it is 
completely decoupled from the medium and thus has no Z dependence. To view in the laboratory 
space-time variables, we recall the definitions of Z and T defined in Eq. ([5]) and write Eq. (146ft as 



az^ =A c ^ + dtr h \-F L )=°- (47) 

Viewed in these variables we see that Poyting's Theorem is still obviously satisfied and that the 
Stokes pulse is moving at the speed of light, due to the fact that it is decoupled from the medium. 

To verify Poynting's Theorem in regime I, we insert the pulse solutions from Eqns. (|36|) and 
the excited state density matrix element from Eq. (|37p into Eq. (|45|) giving 



^ + 2*4<„33> = - ,Z) +^—^±^(1 - K Z). (48) 



Using the definition of k given in Eq. (|28p along with the Doppler averaging function allows us to 
re-write the right hand side of Eq. (|48jl as 

4/1 9 d \ , 2 (T 



thus clearly satisfying Eq. (|4oj) . Unlike the output Stokes pulse the input pump pulse does 
depend on Z. In Poynting's Theorem this dependence is compensated by the corresponding time 
dependence of the excited state density matrix element. Thus when written in laboratory-frame 
coordinates one can see that the pump pulse group velocity is reduced relative to c, due to the 
coherent excitation and de-excitation of the medium caused by the pulse. 

One can similarly verify that the general solutions, given in Eqns. (|26p for the pulses and Eq. 
(|29c|) for the excited state density matrix element, also satisfy Poynting's Theorem given in Eq. 
(|44p for all Z and T and without pure-state assumptions. This shows explicitly that no excitation 
is lost in the pulse transfer through the Raman exchange regime II, despite the initial "obvious" 
implication that it was. 

VII. NUMERICAL RAMAN SOLUTIONS 

Our analytic solutions are clearly highly specialized. To test their broad relevance, we examine 
more "normal" pulse evolution now by numerical methods and look for correspondences with 
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prominent features of the analytic solutions. Different medium preparations and input pulses are 
examined. We use gaussian instead of sech pulses and replace the infinite uniform medium by 
a medium with definite entry and exit faces. In Fig. H] we show three plots, each containing a 
sequence of six snapshots of the pulse evolution. Each plot corresponds to a different medium 
preparation. Given a medium prepared with a 2 > P 2 , the pump must eventually be absorbed and 
the probe pulse amplified, until the pump is fully depleted. It is useful to define Zt to be the 
"transfer length," which is the length of the medium needed for the two pulse Areas to be equal. 
We can estimate this from Eqn. (|4ip by defining 9^—Zt) = 9^ where 9^ is the input area of 
the Stokes pulse (see Fig. [2]). Using this definition, we solve for Zt giving 

* = 53S^ ta ((ipr) (50) 

For the left plot the predicted transfer location is kZx ~ 6, and we see from frame 4 of the left 
plot that this is confirmed. For the center plot we predict nZj* ~ 10, which is also confirmed. For 
the right plot the analytic formula gives kZt ~ 30, but the last frame only shows the pulses to 
kZ ~ 25 so the Stokes pulse is still weak. We note that these confirmations are obtained despite 
the fact that the input pulses are gaussian: 

and the input pump-pulse Area is not 2tt. Because we are starting with a pulse area 9 a < 2tt, the 
actual transfer location is slightly after the predicted location, which also explains the weak nature 
of the Stokes pulse in the right plot of Fig. SJ The most efficient transfer clearly occurs when the 
Raman inversion, a 2 — (3 2 , is greatest. 

The I-II-III three-regime behavior we have described is apparent in all three plots in Fig. 01 
and we find the same behavior for a wide variety of input pump pulse Areas with gaussian shaped 
envelopes and different mixed-state medium preparations. Our analytic solution for the output 
Stokes pulse in regime III, given by Eq. (|34p . is a very good approximation to the numerical Stokes 
pulse. The output pulse tails show the exponential behavior characteristic of a sech shaped pulse, 
and the Area of the output Stokes pulse is nearly independent of the pump pulse Area and near 
to 2tt. That is, even when two pulses copropagate, the single-pulse principles of SIT remain the 
strongest determiner of pulse evolution. 

This conclusion is "universal" in the same sense that 2tt sech pulses are the "universal" conse- 
quences of pulse propagation in two-level media if the input Area satisfies tt < 9 a < 3n (of course, 
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FIG. 4: Snapshots of the numerical pulse solutions of the full three-level Eqs. Q and © for A = 10/r = 
3/T 2 *, and initial pump-probe pulses with Areas 9 a = and Of, — 0.0057T. The horizontal axis is x in 
units of k/c, and the vertical axis is the magnitude of the Rabi frequency in units of r _1 . The shaded zone 
indicates the location of the medium. The solid curve is |fi a | and the dashed curve is |f2f,|. The top, middle, 
and bottom plots differ only in the way the medium is prepared. In all cases a 2 + 1 = 1, and in the three 
plots we have a 2 — /3 2 = 1.0, 0.6, and 0.2, respectively from top to bottom. The results show the emergence 
of an amplified probe pulse that approaches sech shape and Area « 2n, as in SIT. 
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under the fully coherent conditions that the pulse durations are short enough to ignore homoge- 
neous relaxation and the medium is inhomogeneously broadened). For input Area greater than 
this range a more complicated pulse-breakup behavior enters the picture. To illustrate all of this 
concretely we show in Fig. [5] just the input and output pulse shapes starting from input pulses 
that are gaussian with Areas 9 a indicated on the plots and 8), = 0.0057T. In contrast to all previous 
plots, here we plot the pulse envelopes as a function of the retarded time variable T. In the bottom 
frame one can see pulse breakup because the input Area is greater than 3tt. In previous work with 
two-pulse evolution in three-level media [rJ similar breakup was reported in numerical solutions of 
the reduced two-level equations. Here we see that the underlying cause for the breakup is related to 
higher order SIT-soliton behavior in the fundamental three-level equations. We will not pursue the 
study of larger- Area pulses here, but higher order analytic solutions to the fundamental equations 
are available for the case 



VIII. COMPARISON TO ADIABATIC SRS THEORY 



Our short-pulse theory of SRS retains the full three-level character of the medium, but SRS 
is typically studied in the far off-resonant regime where the average detuning A is so large that 
the excited state (level 3) can be adiabatically eliminated from the density matrix equations. The 
remaining two-level adiabatic MB equations have been extensively analyzed in the context of tradi- 
tional Raman scattering with cw pump and probe fields 13] • Both "steady state" (homogeneous 
decay dominant) and "transient" (homogeneous decay neglected) analytic and numerical solutions 
have been found, assuming cw fields. 

Given the analytic solutions presented above, we can examine the adiabatic two-level equations 
in a new way, simply by comparing exact with adiabatic in the limit of large detuning. To quickly 
review the standard adiabatic elimination technique 



35J, HQ, we formally integrate pis'- 



i 

pi 3 = iAp 13 + - [Q a (p 33 - pn) - MbPu], (52) 
with the initial condition pi 3 = at t = — oo, to obtain 

pis(t) = i I e iA (*-*') A(t')dt' '. (53) 



Here A(t) = §[^ a (*)G°33(*) — Pn(*)) ~~ &>b(t)pl2(t)], an d we have suppressed the Z dependences. 
We integrate Eq. (|53|) by parts repeatedly to obtain the usual series 

A(t) i dAit) 1 d 2 A(t) 
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FIG. 5: Plots of numerical pulse solutions of Eqs. (J3j> and ([6]) for "strong" and "weak" input Areas. Here 
the display is given as a function of T instead of Z . The pump pulses have "strong" input Areas as indicated 
on the plots, and the probe pulses have "weak" Areas Ob = 0.0057T in all but the 4th plot. The top frame 
corresponds exactly to the input and output solutions shown in the top plot of Fig. |H The left column of 
snapshots shows the gaussian input pulses at the entrance face to the medium. The right column shows the 
output pulses at the output face after the pulses have propagated nZ = 40 where the output Stokes pulse is 
well described by the sech-shaped 27r output pulse of Eq. (j34|) in all but the last frame. In all frames both 
pump and Stokes pulses are plotted, however due to scale (or overlap as in the 4th plot) just the pump or 
Stokes can be visualized. The last frame shows an example of an SIT-type breakup. All of these solutions 
are for media with a 2 —l and /3 2 =0. 

In the case at hand, where each field consists of a single pulse, the "strength" of the interactions 
can be estimated from the pulse Areas 9 a and 9 b . For the sake of simplest estimates we can take 
both pulses to be roughly similar to the extent that we can denote both of their durations as r, as 
we already indicated in Eqn. ([5T]) . for example. Then we have Q a ~ 9 a /r and f2{> ~ and by 
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extension d£l a /dt ~ 9 a /r 2 , etc. As a consequence 



M , dA M 

~dt 



A(t) ~ — , and ~ — , etc., (55) 



where 



M = ^[9 a (p 33 - Pll ) -9 bPl2 }, 

and M is clearly bounded by the magnitudes of 9 a and 9^, which will be taken of order 2tt or much 
smaller. 

The consequence of these estimates is to turn Eqn. (|54p into a power series in the parameter 
1/rA. In traditional SRS this is a very small parameter because the detuning is typically large 
and the pulse duration is also large or even extremely large in the case of quasi- cw fields. Then 
only the first term is needed in (|54p . and one obtains the familiar "adiabatic following" solution 
/9i3 « — A(t)/A, equivalent to setting P13 = in the density matrix equations. This lowest order 
solution for /J13 and a similar one for ^23 can be paired with the relation ^33 ~ 0, which also follows 
from large rA values when Areas are not large. Then the three-level density matrix equations 
simplify to equations for the density matrix of the two ground states alone: 

Pn = i—^-P2\ - i—^—pn (56a) 

P22 = i—^—pi2 ~ i^-P2i (56b) 

P12 = i—^-[P22 ~ Pu) + «A l V12, (56c) 

where 9.^ = Q, a Q,* h /2lS. is the two-photon Rabi frequency and A^ 2 ) = (\Q a \ 2 - |O b | 2 ) /4A. is the 
two-photon AC Stark shift. The same procedure simplifies Maxwell's equation ([6]), giving 

= -^(^Pii + 6 pi 2 ) (57a) 
= -i^{^bP22 + O a p2i)- (57b) 

In comparing with the analytic solutions we first expand them in a series of powers of the 
small parameter 1/rA, and then insert the lowest order contributions to those solutions into the 
two-level equations. One finds agreement, confirming the lowest-order adiabatic approach, but 
with a significant qualification. An even-handed comparison requires that the two-level equations 
themselves be first reduced to the internally consistent lowest adiabatic order. In the most common 
context, the SRS process involves amplification of a weak Stokes pulse by a strong pump pulse, 
in which case only the leading term in the first of Eqs. (|57p remains, with p\\ = a 2 = constant. 
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This reduces the equation to one that predicts only dispersive (phase) change with no absorptive 
effect at all - no change in pump amplitude. It is straightforward to see that this is consistent with 
the exact solutions to the same order because the propagation coefficient in Eq. ([2"8]) reduces to 
k ~ in the rA 3> 1 limit, preventing any susubstantial changes in amplitude during propagation, 
allowing only phase change. 

The remaining question, therefore, is whether the SRS two-level equations can be confirmed 
beyond this lowest order, the order in which they were derived. What is found by comparison 
beyond lowest order is that the traditional two-level SRS theory does not consistently predict the 
short-pulse Raman transfer process developed in the preceding Section. Important details of those 
solutions are missing. This conclusion can be examined in numerical solutions, as in Figs. [6] and [71 
We first consider injecting two weak pulses, to isolate the effects of a short pulse duration, as shown 
in Fig. [6j The differences in the pulse evolution are clear, showing that short pulse propagation 
is not reliably reproduced. Specifically the SRS two-level equations do not reproduce the reduced 
group velocity and absorption of the pump pulse. Only for short medium lengths is the group 
velocity dispersion negligible, making the SRS two-level equations a valid approximation. 

Further comparison of the SRS two-level equations with the three-level equations can be made 
by examining details of the outputs obtained when we inject two sech pulses with different Areas, 
one of them large and the other small: 9 a = 2ir and 8^ = 0.0057T, and neglect inhomogeneous 
broadening. In the left frames of Fig. [7J we show the pulse envelopes before entering the medium 
and in the right frames we plot the output pulses after exiting the medium at kZ = 40. The first 
row is the solution to the full three- level equations, and the second row is the solution using the 
reduced equations, and the bottom figure is a zoom plot of the output pulses in the second row. In 
the second row one can clearly see the inadequacy of the conventional adiabatic elimination process 
in this case. The Stokes pulse only grows where the pump pulse has been depleted to satisfy the 
Manley-Rowe relation ^(|O a | 2 + \^b\ 2 ) = [easily established via Eqns. ([57]) ]. and the pump pulse 
is never fully depleted. In addition because k ~ in the adiabatic limit, the pulse transfer occurs 
at velocity c. In contrast, the three- level solution accurately fits the analytic transfer solution Eq. 

(USD. 

Incidentally, we note that the complicated pulse shape of the two-level solution in Fig. [7J is 
similar to the early numerical results presented by Tanno, et al. [18], where the reduced two-level 
equations were solved for gaussian shaped input pulses. As a result of our comparisons we can 
now suggest why those earlier numerical solutions resisted clear interpretation. It appears that 
the pulse durations and strengths were used in combinations to make pulse Areas greater than 
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FIG. 6: Snapshots of weak pulse propagation, with horizontal and vertical axis the same as similar plots. Top six 
frames show plots of numerical solution to full three level equations. Bottom six frames show numerical solution to 
reduced two-level Raman equations. The medium in both cases is a mixed state with a 2 = 0.8 and /3 2 = 0.2. Both 
input pulses are weak and sech-shaped with input area 0.005-7T. The three-level equations show both pulses being 
absorbed by the medium with differing group velocities and minimal pump-Stokes transfer as expected for weak 
pulses. The reduced two- level equations cannot accurately describe the reduced group velocity and absorption effects 
that clearly occur in the full three-level case. In this weak pulse example, the two-level equations cause minimal 
amplification of the Stokes pulse, causing the pump and Stokes pulses to appear on the plot as only a single pulse. 



permitted for validity of the twodevel equations. In particular, we believe that the upper level 
should have played a role and needed to be included. 
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FIG. 7: Evolution predicted by the complete three-level equations (top row) and by adiabatic two-level 
equations (second row). The left and right frames show the pulses at input and output before and after the 
medium (at kZ — 40) prepared with a 2 — 1 and f3 2 = 0. The horizontal axis is the retarded time T in units 
of the pulse width r, and the vertical axis measures the product f2r for each pulse. The solid curve is |fi„| 
and the dashed curve is which is hardly visible in the input plots, but is plotted. 

IX. CONCLUSIONS 



We have presented an analysis of the coupled Maxwell-Bloch equations governing two-pulse 
propagation in lambda-type media, under conditions that can be called "short-pulse EIT" . By this 
term we mean that the equations were solved under conditions familiar in EIT scenarios except 
that the pulses are taken very short compared to medium relaxation times instead of very long, 
and both pulses are propagated without typical EIT pump-probe intensity relationships. This is a 
regime where experimental tests have yet to be undertaken, but no fundamental barriers to tests 
appear to be present. 

Our solution formulas for two-pulsejjropagation were obtained by the Backlund transformation 
method as refined by Park and Shin [17], and were built on initial conditions that included inco- 
herently prepared ground states. The solutions are complicated but they have easily understood 
asymptotic forms for both x — ► —oo and x — > +oo, although these limit forms are different from 
each other. This difference is compatible with stimulated Raman scattering, where initial and 
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final stages are very different. The solutions allowed us to undertake comparisons not previously 
made, between propagation in pure-state and mixed-state media, and we could assess the role of 
the atomic "dark state" [20I in a new context. Our results show that there is essentially no role 
here for the dark state as x — > 00, in direct contrast with the results for the same pulse preparation 
in pure-state media fli] ]. 

Our analytic solutions were shown to be stable in the important sense that the output results 
they predict are also obtained from numerical integration of the same equations, but without 
needing to insist on either sech shape or soliton-perfect conditions on the input. This again shows 
that the 2n condition derived in the McCall-Hahn Area Theorem [j| is sufficient to determine 
asymptotic pulse formation in an absorbing medium. We see that the same medium, prepared in 
the same way throughout [—00 < x < +00], serves as a "pump" attenuator in the early stages and 
as a "probe" attenuator in the final stages. 

We took advantage of the numerical approach to make a test of the familiar reduced two-level 
"adiabatic" Raman equations, and found they are not adequate for the short-pulse domain studied 
here. The next term beyond the first-order approximation in the standard adiabatic series must be 
included for agreement with the analytic solution formulas. Of course, in going outside the 2tt Area 
domain one encounters soliton breakup, and this is proposed as the explanation for non-smooth 
pulse shapes obtained in previous numerical studies. 
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